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METHOD AND APPARATUS FOR ESTIMATING 
DOPPLER SPREAD 

FIELD OF THE INVENTION 

The present invention relates to mobile radio communications, and more 
particularly, to estimating Doppler spread in a mobile radio communications system. 

BACKGROUND AND SUMMARY OF THE INVENTION 

In mobile radio communications systems, the signal propagation between a 
transmitter and a receiver can be understood by introducing the concept of a mobile radio 
channel imptilse response that introduces a filtering action on the signal. For example, the 
baseband impulse response, h{t) , can be expressed as: 

where Xi is the amplitude of the i"^ received signal, is the phase shift of the V-^ received 
signal^ 6 is the impulse function, and is the time delay of the i**^ received signal. 
Equation (1) shows that the received signal can be thought of as a series of time-delayed, 
phase-shifted, and attenuated versions of the transmitted signal. If the channel is time 
variant, then A' ^£ ^^^^ functions of time. 

Radio propagation is characterized by multiple reflections, diffractions, and 
attenuations of the signal caused by natural obstacles like buildings, hiUs, etc. which result 
in multi-path propagation. Accordingly, radio waves arrive at a mobile receiver from 
different directions and with different time delays, and are made up of a possible direct ray 
and reflected rays. Scattered reflections cause Rayleigh fading and have very small mutual 
delay differences that are not resolvable in the receiver. These rays (which have a 
magnitude and a phase) combine vectoriaUy at the receiver to give the resultant signal, 
which depends on the phases and amplitudes of the reflections that exist in the multi-path 
field. 
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As a mobile radio moves from one location to another, a Doppler effect 
occurs. In general, the Doppler effect relates to a change in the apparent frequency of a 
source of electromagnetic radiation when there is a relative motion between the source and 
the observer. For moving mobile radios, where multi-path propagation is common, the 

5 phase relationship between the components of the various incoming rays changes, so the 
resultant signal is spectrally spread. This phenomenon is referred to hereafter as Doppler 
spreading. Whenever relative motion exists in the absence of multipath propagation, there 
is a Doppler shift of the frequency components of the received signal, but no spreading of 
the received signal spectrum. This phenomenon is referred to as Doppler shifting. This 

10 present invention is directed to estimating Doppler spread. 

A multi-path signal envelope is charaaerized by a distribution function of 
amplitude that follows the so-called Rayleigh distribution function, which is why mxjltipath 
sometimes referred to as Rayleigh fading. When a mobile radio moves in a multipath 
environment, the received signal appears to vanish, i.e., "fade," at certain positions. 
15 However, moving a few meters brings it back again. Thus, a mobile radio moving in a 
multipath environment experiences signal fluctuations, and this effect is referred to as 
Rayleigh or fast fading. 

Fig. 1 illustrates a schematic diagram of a communications system 10 with a 
transmitter 12, transmitting information over a carrier frequency on a radio channel that 
20 is subject to fading a{t) and noise ?^t). The received signal is received on a frequency that 
is offset from the transmit carrier frequency. The time-varying, complex fading coefficient 
a{t} models fast fading. Li essence, the Doppler spread describes how fast the channel is 
changing, or equivalently the spreading of the received signal caused by Raykigh fading. 
More formally, 



25 



f-L-ic (1) 

where is the Doppler spread, is the carrier frequency, v is the mobile velocity, and 
c is the speed of light. 
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Third generation wireless cellular communications systems, such as the 
Universal Mobile Communications Systems (UMTS), must support commvinications with 
a mobile station traveling at a considerable velocity, e.g., up to 500 km/h, over Rayleigh 
fading radio charmels. The resulting Doppler spreading on such a channel leads to 
5 performance degradations absent suitable compensation. To be able to track the channel 
accurately, knowledge of the Doppler spread is needed. One variable that affects the 
accuracy of such an estimation is a frequency offset caused by differences between local 
oscillators in the transmitter and receiver. 

Fig. 2 illustrates Doppler spreading and frequency offset errors. More 
10 particularly, Fig. 2(a) shows a transmitted signal corresponding to a narrow frequency 
spectrum. Fig. 2(b) shows that transmitted signal spectrum much wider after Doppler 
spreading by a fast fading chaimel. In addition to the Doppler spreading in the frequency 
domain, where shows the maximum Doppler spread of the spectrum, there is a 
frequency error or offset that shifts the whole frequency spectrum. Fig. 3 illustrates phase 
15 changes of the received signal caused by Doppler spreading over a fast fading charmel with 
a frequency offset. The slope (dashed line) of the phase is greater than zero represents a 
positive frequency offset. The impredictable variations in the phase are caused by Rayleigh 
Sading. When a Rayleigh fading signal is received, it can be viewed as a sirni of a number 
of incoming waves, each with its own amplitude, frequency, and phase. The vector stim of 
20 these waves looks in the frequency domain like Fig. 2(b), and in the time domain like 
Fig. 3. In order to accurately determine the Doppler spreading of a fast fading charmel, 
any frequency offset should be compensated before estimating the Doppler spreading. 

Accurate Doppler spread estimation can be quite useful in mobile 
communications systems. First, knowledge of the Doppler spread can be med to improve 
25 the performance of the demodulator and reduce bit error rates. Second, accurate Doppler 
spread estimations may be used to optimize interleaving links in order to reduce reception 
delays. Third, in a cellular system with hierarchical cell layers, the Doppler spread may also 
be used by the network in cell layer assignment strategies. For example, low speed mobile 
radios would be assigned to pico cells, mediiam speed mobile stations to micro cells, and 



high speed mobile stations to macro cells. By assigning high speed mobiles to large cells, 
the number of handovers can be reduced. This reduces the amount of signaling, and 
therefore improves system capacity. 

A fourth example use for accurate Doppler estimation is to adjust the 
filtering bandwidth for charmel estimation in the receiver, which is one way id improve 
demodulator performance. When a narrowband signal is transmitted over a Raylo^ 
fading channel, its spectrum spreads. After this spectral spreading, it looks like Fig. 3 in 
the time domain, hi order to coherently detect the data symbols, a coherent receiver has 
to track these unpredictable changes in the phase. This tracking (often called channel 
estimation) is usually performed using some kind of filtering, with the goal to filter out as 
much noise as possible. With knowledge of the Doppler spread, the filtering can be 
adjusted to pass the useful signal and filter out the noise outside the Doppler spectrum. 

A fifth example use would be to set an appropriate power control step in a 
wideband CDMA-type mobile communications system. For a Rayleigh fading channel 
with a large Doppler spread, the charmel power changes rapidly. Power control tries to 
compensate for these channel changes. Therefore, for a rapidly changbg channel (i.e., 
high Doppler spread), it can be advantageous to use a larger power control step size than 
for a slowly fading channel. 

The present invention achieves these and other benefits by estimating a 
Doppler spread associated with a Rayleigh or fast fading channel established between a 
transmitter and receiver, e.g., a base station and a mobile station. In particular, the 
Doppler spread is estimated through calculation of the autocorrelation function of a 
sequence of complex charmel estimates determined from the known sequence in a 
received signal. More specifically, a sequence of complex channel estimates obtained from 
the known sequence in a first sampling interval is complex conjugated and then correlated 
with a sequence of complex channel estimates obtained from the known sequence in a 
second sampling interval which have been not been complex conjugated. A zero crossing 
of the complex autocorrelation function is detected, and the estimated Doppler spread is 
calculated using this zero crossing and a Bessel function. 
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In a preferred, non-limiting embodiment, the known sequence is 
compensated for a frequency offset. The estimated Doppler spread and the compensated 
known sequence may be used to estimate the Rayleigh or fast fading channel. The 
estimated channel response is used to filter an unknown sequence in the received signal to 
5 compensate for phase error caused by the Doppler spreading over the fast fading charmel. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features, and advantages of the invention 
will be apparent from the following description of preferred, non-limiting example 
embodiments, as well as illustrated in the accompan3^g drawings. The drawings are not 
10 to scale, emphasis instead being placed upon illustrating the principles of the invention. 

Fig. 1 is a function block diagram illustrating a communications system 
where information is transmitted between a transmitter and a receiver over a Rayleigh 
channel; 

Figs. 2(a) and 2(b) illustrate a transmitted signal at the output of the 
15 transmitter in Fig. 1, and the received signal after spreading by the Rayleigh channel, 
respectively; 

Fig. 3 illustrates in the complex plane an example of the received signal 
influenced by frequency offset and Rayleigh channel spreading; 

Fig. 4 is a block diagram of a mobile radio communications system which is 
20 one example context in which the present invention may be advantageously employed; 

Fig. 5 illustrates a format of transmitted data which includes both known 
and iinknown sequences per timeslot; 

Fig. 6 illustrates a portion of a receiver to illustrate using analog components 
how a received complex signal may be complex conjugated; 

25 Fig. 7 illustrates an example Doppler spreading estimator; 
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Fig. 8 illustrates one example of autocorrelation in accordance with the 
present invention; 

Fig. 9 illustrates an autocorrelation function in accordance with the present 
invention as a function of lag identifying a first zero crossing of that autocorrelation 
5 function; 

Fig. 10 illustrates an interpolation technique for determining that zero 

crossing; 

Fig. 11 illustrates a block diagram of an example receiver in which the 
;l present invention is employed; and 

10 Fig. 12 illustrates example procedures for estimating Doppler spread. 

\M 

i " DETAILED DESCRIPTION OF THE DRAWINGS 

i . In the following description, for purposes of explanation and not limitation, 

I' specific details are set forth, such as particular embodiments, procedures, techniques, etc., 

; in order to provide a thorough understanding of the present invention. However, it will be 
~ 15 apparent to one skilled in the art that the present invention may be practiced in other 
embodiments that depart from these specific details. In some instances, detailed 
descriptions of well-known methods, interfaces, devices and signaling techniques are 
omitted so as not to obscure the description of the present invention with tmnecessary 
detail. Moreover, individual function blocks are shown in some of the figures. Those 
20 skilled in the art will appreciate that the funaions may be implemented using individual 
hardware circuits, using software funaioning in conjunction with a suitably programmed 
digital microprocessor or general purpose computer, using an Application Specific 
Integrated Circuit (ASIC), and/or using a Digital Signal Processor (DSP). 



The present invention may be employed in any mobile radio 
communications system such as that shown in Fig. 1. However, a more specific, but non- 
limiting example of a mobile communications system is the Universal Mobile 
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Telecommunications System (UMTS) 10 shown in Fig. 4. An example Universal Mobile 
Telecommunications System (UMTS) 10 is shown in Fig. 1. A representative, core 
network 16, includes a circuit-switched core network (CS CN), shown as box 18, and a 
packet-switched core network (PS CN), shown as box 20. The circuit-switched core 

5 network includes nodes, e.g., Mobile Switching Centers (MSC) 18, Home Locations 

Register (HLR), Gateway MISC (GMSC), etc., that provide circuit-switched services. The 
packet-switched core network includes nodes, e.g., Serving GPRS Support Nodes 
(SGSN) 20, Gateway GPRS Support Node (GGSN), HLR, etc., that are tailored to 
provide packet-switched services. The CSCN 18 is coupled to an external circuit-switched 

10 network 12, such as the Public Snatched Telephone Network (PSTN) or the Litegrated 
Services Digital Network (ISDN). The packet-switched core network 20 is coupled to an 
external packet-switched network 14, such as the Internet. 

Each of the core networks 18 and 20 is coupled to a UMTS Terrestrial 
Radio Access Network (UTRAN) 22 that includes one or more Radio Network 

15 Controllers (RNCs) 26. Each RNC is coupled to a pliu-ality of base stations (BSs) 28 and 
to other RNCs in the UTRAN 22. Each base station 28 corresponds to one access point 
(one sector or cell) or includes plural access points. Radio commtinications between one 
or more base station access points and wireless user equipment unit (UE) 30 are byway of 
a radio interface. Radio access in this non-limiting example is based on Wideband-CDMA 

20 (W-CDMA) with individual radio channels distinguished using spreading codes. 

Wideband-CDMA provides wide radio bandwidth for multi-media services including 
packet data applications that have high data rate/bandwidth requirements. While Doppler 
spread estimation may be employed in any radio receiver, including the UE 30, it is 
particularly advantageously employed in radio base stations 28. 

25 The Doppler spread estimation technique in accordance with the present 

invention may be used for both uplink (mobile-to-base station) and downlink (base 
station-to-mobile) communications. The example slot structure described herein is similar 
to a slot structure used in the uplink in W-CDMA In the downlink, a common channel, 
such as the Common Pilot Channel, may be used for channel and Doppler estimation. 
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Referring to Fig. 5, radio transmitters, such as user equipment units 30, 
typically transmit information over a Rayleigh or fast fading channel to base stations 28 in 
timeslots indicated as . From the known sequence in each time slot, a complex 
sequence of chaimel samples is generated, where n is the time slot index. The unknown 

5 sequence which pertains to the actual content to be communicated to the user equipment 
unit 30 is sent to a demodulator for detection. The present invention employs the known 
sequence in each timeslot to estimate the Doppler spread associated with the Rayleigh or 
fast fading charmel between the user equipment 30 and the base station 28. Since every 
user equipment unit moves at a user determined velocity, each channel has its own 

10 Doppler spread. As explained above, Doppler spread is defined as a spreading of the 
transmitted signal spectnom in the presence of multipath propagation and movement of 
the user equipment unit. Doppler spread is a measure of the rate of the random phase 
variations introduced by the Rayleigh channel. 

While Doppler spread can be estimated in different ways, the present 
15 invention estimates Doppler spread by calculating a complex autocorrelation funaion for 
a sequence of complex channel estimates. These complex channel estimates are phasors or 
vectors in the complex plane with a magnitude and a phase angle (as opposed to a scalar 
which has only a magnitude and no phase angle). The Doppler spread estimation is based 
upon the principle that the autocorrelation funaion of a Rayleigh or fast fading channel 
20 Q may be expressed as a zero order Bessel function of the first kind, /q as set forth 
below in equation (3): 

Q=/o(2^-/.-/-7;) (3) 
where / is the "lag" in the autocorrelation function Q . The lag in an autocorrelation is 
the time shift or difference between correlated signals. More specifically, the lag is the 
25 mutual time difference when two versions of a sequence are correlated. For example, the 
autocorrelation for lag = 0 for the sequence x would be 



C(0) = X (O) X (O) + X (l) X (l) + X (2) X (2) + . . 
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and for kg = 1 



C(l) = X (0) X (1) + X (1) X (2) + X (2) X (3) + .. . 



(in tliis example x is a real signal). 
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At a particular lag value, , where the autocorrelation Q is zero, 
equation (3) may be set to zero and rewritten to solve for the Doppler spread as foUows: 



271 -i^ -T^ 



(4) 



where is the estimated Doppler spread in Hertz, is the lowest positive value that 



where /□ is the zero-order Bessel function of the first kind. 

Using these relationships, the present invention determines the 
autocorrelation function of a received signal using the known sequences, detects the first 
zero crossing of that autocorrelation function, and solves for the Doppler spread using the 
15 determined zero crossing lag 4 and equations (4) and (5) above. 

The autocorrelation for a received signal is determined based on several 
timeslots of known sequences (see Fig. 5). The number of timeslots used for estimation is 
a design parameter that depends on the desired accuracy of the estimate. The more 
samples used, the more accurate the estimation, but at a cost of increased computational 
20 complexity. A typical value for this is 60, but it depends on input signal to noise ratio. 
The input to the Doppler spread estimator in accordance with the present invention is a 
sequence of complex channel estimates {p„ ) acquired from the known sequences 
transmitted. The autocorrelation function of this known sequence is determined in 
accordance with: 



satisfies equation (5) below, 7] is the sample time interval corresponding to one sample 
per timeslot, and 
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(5) 
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c, = j:ip.)-(p..r (6) 

k=0 

where / is the lag in the autocorrelation function, and N is the number of channel 
coefficients used for estimation. An asterisk (''") denotes the complex conjugate. The 
autocorrelation in the present invention uses complex values so that phase changes as well 
5 as magnitude changes are taken into account when determining an estimate of Doppler 
spreading. 

Fig. 6 illustrates a simplified portion of a radio receiver showing one 
example of how complex conjugates of complex sequences may be obtained using analog 
components. The received radio signal (which is real) is obtained by way of antenna 40 

10 and amplifier 42 before being provided to two mixers 44 and 46. The other input to 

mixer 44 is the cosine of the carrier frequency generated by the local oscillator 48. The 
other input to mixer 46 is the sine (90° out of phase with the cosine) of the carrier 
frequency generated by the local oscillator 48. The output of mixer 44 is the real, in- 
phase portion (/) of the sequence, and the output of the mixer 46 is the imaginary or 

15 quadrature portion (Q) of the complex sequence. The imaginary portion is multiplied by 
the imaginary unit "/' represented here as amplifier 50. The in-phase and quadrature 
components are added in a summer 54 to obtain the complex number / + jQ. Complex 
conjugation is more commonly and more preferably performed in the digital domain after 
analog-to-digital conversion, because it is a simple operation and does not require the 

20 expensive hardware depiaed in Fig. 6. The analog implementation shown in Fig. 6 is only 
for purposes of illustration. 

To obtain the complex conjugate of a complex number, the sign of the 
imaginaiy nxamber is changed, i.e., / - jQ. This change in sign is represented by the 
inverter 52. Thus, when a complex conjugate is desired, the inverter path is selected; 
25 otherwise, the quadrature sequence is passed without inversion to summer 54. The 

complex number output from summer 54 is then passed to a sampler or analog-to-digital 
converter 56 which samples the complex sequence once every sampling interval 7| . Since 
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some of the symbols in the sequence are known to the receiver, those known symbols can 
be multiplied with the samples of the received signal corresponding to the known signal in 
order to estimate the effect of the channel. 

Fig. 7 shows one example of a Doppler spread estimator in accordance with 
5 the present invention. The input to the complex autocorrelator 60 is complex channel 
samples generated by a sampler/A-to-D converter 56 shown in Fig. 7. Each of the black 
dots shown as an input corresponds to a sequence of known symbols representing one 
complex channel estimate. The complex autocorrelator 60 calculates the autocorrelation 
function in equation (6). Fig. 8 illustrates a way of performing the complex correlation of 
10 equation (6). The complex autocorrelator 60 uses one charmel estimate from every 

timeslot , and in this example, performs a complex correlation with the known signal for 
a nximber of different lags. For example, the complex correlation for lag 1 is 

=Po- Pi* +Pr Pi* +P2-P3* ■ + ■■■■ 
The complex correlation for lag 2 is 

15 C^=Po-P2*+Pl-P3*+P2-P4*+-- 

The output of the complex autocorrelator 60 is the autocorrelation function C, illustrated 
in Fig. 9. 

The zero crossing detector 62 determines the first zero crossing of the 
autocorrelation function. At this stage, the estimated autocorrelation function should be 

20 compared with the theoretical one, i.e., the Bessel function. This can be done by finding 
the location of the first zero on the positive axis and comparing it to the first zero in the 
estimated autocorrelation function. Locations of the first minimum in the autocorrelation 
function could also be compared. The search for the first zero crossing is performed by 
identifying the first lag value for which the autocorrelation function is negative. Fig. 9 

25 illustrates particular lag values for the autocorrelation function and indicates the first zero 
crossing at lag . That zero crossing does not correspond to an integer lag value, and 
therefore, the exact lag value corresponding to the zero crossing may be determined by 
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linear interpolation as shown in Fig. 10. For example. Fig. 10 shows lag 6 and lag 7 as the 
integral lag values adjacent to the zero crossing of the autocorrelation function. Therefore, 
the zero crossing lag value can be linearly interpolated to be 6.5. 

Returning to Fig. 7, the zero crossing lag value is then provided to the 
5 Doppler spread calculator 64 which solves equation (4) above using the value, the 
known timeslot/sampling duration , and equation (5). 



Fig. 11 illustrates an example radio receiver 70 in which a Doppler spreading 
estimator 80 in accordance with the present invention may be advantageously employed. 
A received signal is parsed so that the known sequence is input into a mixer 74, and the 

10 unknown sequence is input to a mixer 86. Before inputting the known sequence p„ to the 
Doppler spreading estimator 80, it is desirable to compensate the known sequence for a 
frequency offset between the local oscillators and the transmitter and the receiver. A 
frequency offset on the signal at the input to the Doppler spreading estimator 80 means 
that its autocorrelation is no longer a Bessel function, thereby impairing the accurate 

15 determination of the Doppler spread. The frequency error estimator 78 determines the 
frequency offset which is amplified at 82 and integrated at 84. The frequency offset is 
employed by a phase rotator 72 which generates a rotating phasor by calculating e^^'^ , 
where / is the frequency of rotation and t is the time. The compensation phasor is 
complex conjugated before being multiplied in mixer 74 with the known sequence . 

20 The frequency compensation performed by mixers 74 and 86) multiplies the phasors from 
the phase rotator 72 with the input sequence. To compensate the input signal rotating 
with the same frequency , the frequency /j is estimated as frequency f2 . Because the 
rotation of the input signal may be expressed by e^^'^'' and the rotation of the 
compensation signal as e''^'^'^ ,the output from the multiplier wiU be 

25 e^^'^^' ■{e'^'^^') = e'^''^^^~^^^\ If is close to /j , the frequency on the input sequence is 
compensated. 
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In any event, the known sequence (preferably frequenq^ offset 
compensated) is provided to the Doppler spread estimator 80 which performs the Doppler 
spreading estimation as described above. The Doppler spread estimate is provided to the 
channel estimator 76 along with the known sequence used to estimate the channel The 

5 channel estimator 76 maybe, for example, a finite impulse response (FIR) transversal filter 
used to estimate the channel impulse response h{t) (see equation (1)) based on the known 
sequence. In accordance with the present invention, that channel estimator takes into 
account the Doppler spread when formulating the channel estimate. Accordingly, when 
the unknown sequence is received, and preferably frequency offset compensated at 

10 mixer 86, it is combined with the channel estimate at mixer 88 to generate a phase- 
compensated sequence. In effect, the unknown sequence is filtered by the updated 
channel estimate. The phase compensated sequence is then forwarded to a demodulator 
to extract the transmitted data content. Depending on system performance requirements, 
the estimated variables m^ be used for different purposes in the demodulator. For 

15 example, a multi-path searcher in a W-CDMA system, used to find new multipaths, may 
use the estimations of Doppler spread and frequency offset to make path searching more 
efficient. 

Fig. 13 illustrates a Doppler spread procedure (block 100) in accordance 
with one example embodiment of the present invention. A transmitted signal having 

20 timeslots including known and unknown sequences is received over a Ra}deigh fading 
channel (block 102). The received signal of the known sequence is sampled (block 104), 
and preferably compensated for any frequency offset (block 106). Complex channel 
estimates are obtained from the known sequence in various time intervals/ slots 
(block 108). An autocorrelation function is calculated using complex channel estimates 

25 and their complex conjugates, see equation (6) (block 110). From the determined 

autocorrelation function, a first zero crossing is detected, e.g., by interpolation (block 112). 
A Doppler spread associated with the Rayleigh fading channel is estimated using the lag 
value for the first autocorrelation function zero crossing and the zero order Bessel 
function of the first kind (block 114). The estimated Doppler spread is used to adjust the 
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channel estimate, and perhaps to aid in other tasks Hke cell assignment, mobile radio 
location, etc. as discussed above (block 116). 

Thus, the present invention accurately estimates Doppler spread, which may 
be used in a number of advantageous ways, some of which are described above. Because 
complex channel estimates are used, the Doppler spread estimation is not sensitive to 
variations in transmit power present, for example, in a power controlled W-CDMA system. 
Since the Doppler spread estimator is placed inside a frequency offset tracking loop, it is 
not sensitive to frequency offsets often present in communication systems. 

While the present invention has been described with respect to particular 
example embodiments, those skilled in the art wiU recognize that the present invention is 
not limited to those specific embodiments described and illustrated herein. Different 
formats, embodiments, adaptations besides those shown and described, as well as many 
modifications, variations and equivalent arrangements may also be used to implement the 
invention. While the present invention is described in relation to a preferred example 
embodiment, it is to be understood that this disclosure is only illustrative and exemplary of 
the present invention. Accordingly, it is intended that the invention be limited only by the 
scope of the claims appended hereto. 



